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Abstract 

A nascent neutron star may be exposed to fallback accretion soon after the proto-neutron star 
stage [l-3]. This high accretion episode can submerge the magnetic field deep in the crust jj- 
The diffusion of the magnetic field back to the surface will take hundreds to millions of years 
depending on the amount of mass accreted and the consequent depth the field is buried [l\ . Neutron 
stars with large kick velocities will accrete less amount of fallback material [8[ leading to shallower 
submergence of their fields and shorter time-scales for the growth of their fields. We obtain the 
relation Tohm oc v ~ 3 / 4 between the space velocity of the neutron star and Ohmic time-scale for the 
growth of the magnetic field. We compare this with the relation between the measured transverse 
velocities, v± and the field growth time-scales, n/fi, inferred from the measured braking indices. 
Although the number of data is small our results clearly support the view that the magnetic field 
of pulsars grow as a consequence of field diffusion to the surface following fallback accretion. 

Keywords: Individual Pulsars; PSR B0531+21 (Crab Pulsar), PSR B0833-45 (Vela Pulsar), PSR B0540-69 , 
PSR J0537-6910, PSR J1833-1034 —stars: neutron —stars: magnetic fields 



INTRODUCTION 



Soon after the discovery of radio pulsars 



9J their nature as rapidly rotating highly magne- 



ized neutron stars (NSs) radiating at the expense of their rotational energy was established 



10j. A dimensionless parameter related to the spin-down torques on these rotationally pow- 



ered pulsars (RPPs) is the braking index defined operationally as n = vv/v 2 , where v is the 
spin frequency and dots represent time derivatives. The value of n should be 3 if RPPs are 
spinning down with magnetic dipole radiation. Most of the measured pulsar braking indices 
liHl4j are close to 3 but slightly less as shown in Table 1 which indicate that another pro- 



cess also contributes to the MDR in braking these objects. The recently measured braking 



index of PSR J1734-3333 as n = 0.9 ± 0.2 



15J and that of J0537-6910 as n = -1.5 [16| 



together with the earlier measurement as n = 1.4 ± 0.2 of the Vela pulsar [17J imply that 
this process might severely alter the spin history of young neutron stars. 

3EJ 



Different mechanisms have been invoked 



indices of RPPs less than 3. An early suggestion 2l| which has been recently advocated 



15 



22] is that the braking indices of RPPs are less than 3 because their magnetic fields are 



growing in time. If the magnetic dipole moment of a pulsar is changing in time, the braking 
index becomes 

n = 3 + 2^- , (1) 

fl v 

where fi is the magnetic dipole moment of the NS. The characteristic time-scale for the 
growth of the magnetic dipole moment, r M = is then 



where r c = —vj1i> is the characteristic spin-down age. 



A RELATION BETWEEN KICK VELOCITY AND FIELD GROWTH TIMESCALE 

A common ingredient of supernova models is the fallback of some material that can not 
reach the escape velocity ffl. Acco rd i„ g to Chevalier the faHbaeh due to the reversed 
shock in SN87A reached to the surface of the neutron star 2 hours after bounce jsj. The 
initial accretion rate of fallback matter can be very large. For SN87A the initial accretion 
rate is estimated ^| to be 35OM yr _1 . The amount of matter that reaches the surface and 



TABLE I. Pulsars with accurately measured braking indices. 



Pulsar 




V 




n 




v± 


References 






(Hz) 


(kyr) 




(kyr) 


(km s ) 




B0531+21(Crab) 


30.225 


1.2399 


2.51(1) 


10.1(2) 


140±8 


[11, 25] 


DUOOO 


^±o y v eid j 


11.2(5) 


11.303 


1.4(2) 


28(4) 


62±2 


[17, 26] 


J1833- 


-1034 


16.159 


4.8535 


1.8569(6) 


16.984(9) 


125±30 


[14, 26] 


B0540- 


-69 


19.738 


1.6763 


2.087(7) 


7.34(6) 


1300±612 


[26, 27] 


J0537- 


-6910 


62.038 


4.9743 


-1.5 


4.422 


634±50 


[16, 26] 


B1509- 


-58 


6.6115 


1.5648 


2.832(3) 


37.3(7) 


* 


[28] 


J1846- 


-0258 


3.0743 


0.72736 


2.65(1) 


8.3(2) 


* 


[12] 


J1119- 


-6127 


2.4512 


1.6078 


2.684(2) 


20.4(1) 


* 


[13] 


J1734- 


3333 


0.85518 


81.280 


0.9(2) 


0.16(2) 


* 


[15] 



Numbers in parenthesis are last digit errors. 



the time required for this to happen varies depending on many details like the initial density 
of the surrounding medium, relative velocity of the NS with respect to the medium, the 
magnetic field and spin frequency of the NS before fallback 23]. 



The rapid accretion of fallback matter can bury the preexisting magnetic field that was 



formed in the proto-NS stage jJ-7, |24|. The field would then diffuse back to the surface in 
an Ohmic time-scale rohm — (A-R) 2 : /f] where AR ~ 0.1 km is the depth of submergence and 
7] = c 2 /4tt<j is the magnetic diffusivity. The ohmic time-scale varies within a large range, 
hundreds to millions of years, depending on the conductivity, a of the crust and AR which, 
in turn, depends on the amount of fallback, AM. If the initial accretion rate is very large 
the magnetic field is submerged rapidly and the amount of accreted matter is determined 
by the density of the medium and the velocity of the NS. 

For a depth at which the density is greater than 10 6 g cm -3 the electrons in the crust are 
ultra- relativistic and the relativity parameter xp = V?l m cC ^> 1 where pp is the Fermi 
momentum and m e is the electron mass. In this regime conductivity set by impurity 
scattering, which is expected to prevail at most of the lifetime of observed pulsars, is 
a = 4.26 x 10 21 x-pZ/Q s^ 1 29[ where Z is the atomic number and Q ~ 1 is the impu- 
rity parameter. To a good approximation the mass AM and thickness AR of the accreted 



layer is a small fraction of the total mass, M and radius R of the star, respectively. Within 
his approximation the relativity parameter depends on z = AR/Hr as xf = \J z(z + 2) 



30(] where i? R = m e c 2 / g p e m p is the scale-height. Here 5- = GMR~ 2 (1 - 2GM/Rc 2 y\ 
fi e = A/Z (A is the atomic weight) and m p is the proton mass. For M = 1.4M Q , i? = 10 
km and p e = 2 one obtains i?R = 7.25 m. The mass of the crust region of thickness AR 
is AM ~ AirR 2 J^_/^ R p( r )dr where p is the density. By using the hydrostatic equation 31] 
this can be written as AM = (4nR 2 / g)P^ where Pb is the pressure at the bottom of the 
layer. For ultra-relativistic electrons P = PqX f where Pq = (27r/3)m e c 2 (m e c/h) 3 = 1.2 x 10 23 



dyne cm 2 . If one approximates \/z(z + 2) ~ z as appropriate for > 1, then 

/ gAM \ 1/4 / AM \ 1/4 

relates the accreted mass to the thickness of the resulting layer. With the same approxima- 
tion rohm oc x F cx (AR) 3 and we find 

/ AM \ 3/4 

r Ghm = 31.5 kyr (Z 26 /Q) [^z^) (4) 

where Z 2 q = Z/26 appropriate for iron. 

The newborn neutron star will be accreting from a uniformly expanding medium where 
density decreases as p oc t~ 3 [3]. The neutron star moving with velocity v in this medium is 
subject to a mass inflow rate of M m = 4tt(GM) 2 pv~ 3 [8(. We combine these together as 

M in = M Q fl + -J , M = -^-^ (5) 

where po is the initial density of the medium and to is the dynamical time-scale at which 
the accretion rate declines. We use the symmetry of the fluid dynamic equations under 
translations in time to get rid of the singularity at t — 
rate can be very high and the amount of matter accreting onto the neutron star will be 
Eddington limited (M* = M E ) until h = t [(M /M E ) 1 ' 3 - 1] at which M m drops down to 
Me = L F R/GM where L F = 4nGMm p c/aT is the Eddington luminosity. 

The accretion rate onto the neutron star, for t > t\ 3> to is M* = M F (t / 1\)^ 3 . This 
implies that the total mass accreted onto the neutron star will be AM = |Mgti. Assuming 
M > M E we find AM = f M e £ (M)/Me) 1/3 oc v' 1 and this leads to 

Tohm = 20 kyr ( lQi2 gr°cm-3 s s ) ( lQOkL-0 ^ (6) 



321 ] . Initially, the mass inflow 



where we scale the poorly known quantities po and to together. 

As a test of this relation rohm oc v _3 ^ 4 we check the relation between inferred from 
braking indices via Eqn.pj and measured transverse velocities v± expecting r M ~ Tohm 
and v± ~ v. Only 9 RPPs have accurately measured braking indices that would allow for 
inferring the time-scale for the growth of their magnetic moment via Equation (J2J). Of 
these RPPs, we could find only 5 with measured transverse velocities as shown in Table 1. 
In Figured] we plot versus v± of these objects together with the estimate of the Ohmic 
time-scale in Eqn.([6]) depending on the space velocity v. Although the number of data is 
small, an inverse relation between the field growth time-scale and the transverse velocities 
can still be noticed as a result of the strong dependence of accretion rate on v. The scattering 
in the data is most likely the result of the weak dependence on the initial density po of the 
surrounding medium and less likely the initial spin and magnetic field of the neutron star 
as these would not play a significant role once the field is buried. Note also that does not 
necessarily coincide with rohm if the field growth is not exponential. Add to this that v± 
may not be a good estimate of the space velocity v if the radial component of the velocity is 
not negligible compared to v±. Given these deteriorating factors the relation obtained with 
5 data can be considered satisfactory. 



DISCUSSION 

As a prediction of the fallback submergence and post-fallback growth of magnetic fields, 
we have derived the relation rohm ck v ~ 3//4 between the space velocity v and Ohmic time- 
scale for the diffusion of a buried magnetic field, Tohm- We found a similar relation between 
the measured transverse velocities v± and characteristic field growth time-scales r M = \xj \i 
inferred from the measured braking indices. Such an inverse relation between the kick veloc- 
ity and the field-growth time scale of pulsars is a prediction of only post-fallback magnetic 
field growth model and not of any other model. 

Our estimations indicate that the depth of field submergence for the considered population 
of young neutron stars is very shallow, possibly leading to the selection of these objects as 
rotationally powered pulsars. Smaller space velocities result with larger amount of accretion 
and very long diffusion time-scales of buried magnetic fields leading to young neutron stars 
not appearing as rotationally powered pulsars. 

5 



Depending on the angular momentum distribution in the progenitor the fallback matter 



can have sufficient angular momentum to form a disk 



33( after the brief high accretion 



episode. Accretion from such a disk is proposed to explain the persistent X-ray luminosity 



of anomalous X-ray pulsars 



34j and all enigmatic young neutron stars [35|. Detections of 



debris disks around magnetars 4U 0142+61 36j and IE 2259+586 [37| motivate further work 
for investigating whether fallback disks are the key ingredient for explaining the diversity 
of neutron stars. As the progenitors of magnetars are at the highest mass and angular 
momentum end of the spectrum their fallback disks might also be the longest lifetime end 
of the spectrum. 
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FIG. 1. The relation between the measured transverse velocities v± and characteristic magnetic 
field growth time-scale = p,/fi inferred from measured braking indices via Equation ([2]) for the 
5 pulsars, B0531+21(Crab), B0833-45(Vela), J1833-1034, B0540-69 and J0537-6910. These 5 
objects form the subset of pulsars with accurately measured braking indices and measured trans- 
verse velocities. Also shown on the plot is the estimate of the Ohmic time-scale given in Eqn.JB]) 
depending on the space velocity v: Dotted line stands for rohm = 358 kyr fu/(km s -1 )) 3 ^ 4 = 
11.3 kyr (v/(100 km s -1 )) 3 ^ 4 corresponding to po i[j = 10 11 gr cm -3 s 3 . Dashed line stands 
for rohm = 636 kyr (w/(km s -1 )) 3 ^ 4 = 20.1 kyr (y /(100 km s -1 )) 3 ^ 4 corresponding to 
Po = 10 12 gr cm -3 s 3 . Dashed-dotted line stands for rohm = 1130 kyr (y/(km s -1 )) 3//4 = 
35.7 kyr (W(100 km s -1 )) 3//4 corresponding to po ijj = 10 13 gr cm -3 s 3 . 
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